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Sunto Viene brevemente descritta la teoria di Black & Scholes che nel 1973 
hanno considerato il problema della determinazione del prezzo di opzioni di 
tipo europeo. Tale procedimento viene successivamente applicato alle opzioni 
asiatiche; in questo caso il problema si formalizza in un problema di Cauchy 
per un’equazione differenziale di tipo parabolico degenere, per il quale vengono 
presentati alcuni risultati di esistenza, unicità ed approssimazione numerica. 


Abstract We briefly recall the basic ideas of the work by Black & Scholes 
that considered the problem of the pricing of an european option. We next 
show that, when we apply the same technique to some asian options, we find 
a Cauchy problem for a degenerate parabolic partial differential equation. We 
present some existence and regularity results for such a Cauchy problem and 
some finite difference schemes for its numerical solution. 
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1 Introduzione 


Presenterò alcuni problemi che hanno origine in Matematica Finanziaria, che mì 
sono stati proposti da Emilio Barucci e Vincenzo Vespri ed hanno dato origine 
ad una collaborazione di cui esporrò i primi risultati. Vorrei anche ricordare che 
le prime questioni di tipo modellistico e probabilistico sono state approfondite 
nel corso della preparazione della Tesi di Laurea di Benedetta Santarelli. 


I problemi che abbiamo considerato consistono nello stabilire un prezzo equo 
per le opzioni, che sono contratti che assicurano al possessore il diritto di ac- 
quistare (o vendere) in futuro un determinato pacchetto di titoli ad un prezzo 
concordato. 

Prima di parlare più dettagliatamente di questo tipo di prodotti finanziari 
vorrei fare alcune osservazioni. Innanzitutto il possesso di un’opzione non corri- 
sponde ad una effettiva ricchezza, ma solamente ad una opportunità futura. No- 
nostante ciò, le opzioni sono titoli attualmente molto trattati: infatti la Chicago 
Broad Options Exchange, pur trattando solo opzioni, è tra le più i mportanti 
borse del mondo ed è seconda solamente alla New York Stock Excange. 

Una prima, fondamentale, risposta al problema della valutazione del prezzo 
di un particolare tipo di opzione è stata fornita nel 1973 da Black e Scholes [4]; 
altri risultati in questa direzione sono dovuti a Merton, con un metodo che si 
è rivelato applicabile ad un’ampia classe di problemi finanziari. Questi studi 
hanno avuto un’enorme ricaduta: da un lato i criteri oggettivi per la valuta- 
zione di titoli rischiosi ha contribuito ad accrescere i mercati, da un altro lato si 
è sviluppato un ampio campo di ricerca di matematica applicata all'economia. 
L'importanza di tale ricaduta è testimoniata dal fatto che nel 1997 è stato asse- 
gnato a Scholes e Merton il premio Nobel per l'Economia. Il merito principale 
di questa teoria è stato tuttavia riconosciuto a Black, deceduto nel 1995. 


Nel seguito di questa introduzione descriverò brevemente l’analisi di Black e 
Scholes, che riguarda le opzioni dette “di tipo europeo”, nel paragrafo successivo 
verranno considerate le opzioni di tipo asiatico, che costituiscono una variante 
di quelle europee. Nei due paragrafi successivi verrano studiati i problemi ma- 
tematici legati a questi ultimi titoli. 


Nel lavoro di Black & Scholes [4] viene considerato il problema della valuta- 
zione di un’opzione europea di tipo “call”. Si tratta di un contratto che assicura 
il diritto, ma non l’obbligo, di acquistare in una certa data T e ad un prezzo 
concordato E un pacchetto di azioni il cui valore verrà indicato con S(t). 

Si assume che S sia descritto dalla seguente equazione differenziale stocastica 


dS = SodW + Sydt (1.1) 


dove o e pi sono fissate costanti positive e W è un processo di Wiener. 

Il valore dell'opzione, che verrà nel seguito indicato con V, è perfettamente 
noto all'istante t = T, quando sarà noto il valore del pacchetto S(T), ed è 
chiaramente 


V(T, S(T)) = max{S(T) — £,0}. (1.2) 
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Black & Scholes hanno proposto un metodo per ricostruire il valore V(t, 5) per 
i tempi precedenti t < T, basato oltre che sulla teoria stocastica, anche sulla 
cosiddetta “condizione di non arbitraggio”, che è un principio secondo il quale il 
mercato non offre agli investitori l’opportunità di ottenere un profitto immediato 
da certe operazioni non rischiose. 

Se si considera V come funzione delle due variabili indipendenti (t, 5) allora, 
per il lemma di Itò, deve risultare 


(1.3) 


ds d6ta°% + 


2 
dV = os awd (st 1 2g2 DV 3) dt. 
Sarebbe opportuno eliminare la componente aleatoria del processo stocastico. 
A tal fine si costruisce un “portafoglio” costituito da una quantità di opzioni, 
dal valore V, e da una quantità A di azioni dal valore S. Indicheremo con 
II=V_-AS il valore di tale portafoglio. Qui si introduce una nuova ipotesi, di 
carattere modellistico: si suppone che A possa variare nel tempo, ma che il suo 


differenziale stocastico dia un contributo trascurabile. Pertanto si avrà 
dIl = dV — AdS, 
ossia, utilizzando le formule (1.1) e (1.3) 
oV OV 1 0°V av 
dl = — — — RT * E 
o$ (3 A)aw+ (us (54 A) 430% 952 + e) di (1.4) 


Se scegliamo 


xa (1.5) 


otteniamo l’equazione “deterministica” 


_ (Laga 884. 0V 
dm= (3o°90t, dt. 


Se ora pensiamo di investire lo stesso capitale in azioni non rischiose otteniamo 
dIl = rIldt 


(qui r indica il tasso di interesse del titolo non rischioso), per il principio di non 
arbitraggio questi due ultimi investimenti devono fornire lo stesso rendimento, 
di conseguenza si ottiene l’equazione di Black & Scholes 
Il aga0V_ 0V OV 
30°S = +rS=-rV=0. 1.6 
o gt gg i. 
In questo modo il problema della valutazione di un’opzione di tipo euro- 
peo è stato trasformato in un’equazione alle derivate parziali di tipo parabolico 
retrogrado (1.6), che appare degenere per la presenza del coefficiente S? che 
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moltiplica la derivata DE, a cui è associata la condizione finale (1.2). Osser- 
viamo che il problema appena descritto è equivalente al più familiare problema 
di Cauchy per l’equazione del calore 


{ Ur = Uro per (2,t) € Rx]0, T{, 


u(2,0) = (2) perzeR, (1.7) 


come si verifica effettunado il cambiamento di variabile x = log(5) e qualche 
aggiustamento secondario e ponendo 


2r- 0? 4r°+40?2r +04 ox 
(x) = exp (il: + frste tor) max {exp (5) _ E,0} ‘ 


Per questa scelta del dato è addirittura possibile scrivere esplicitamente la so- 
luzione di (1.7) per mezzo di una espressione che viene correntemente detta 
i formula di Black & Scholes. 


2 Opzioni di tipo asiatico 


Le opzioni europee costituiscono il più semplice tipo di opzioni trattate; in 
effetti esistono innumerevoli tipi di contratti di questo genere, ognuno dei quali 
dà origine a problemi matematici completamente diversi. 

Ad esempio, le opzioni “americane”, che costituiscono 180 % di tutte le 
opzioni trattate, si caratterizzano semplicemente per il fatto che il possessore 
può esercitare il diritto a vendere (o ad acquistare) il pacchetto in un momento 
qualsiasi t € [0,7] e non necessariamente alla scadenza (t = T) e danno origine 
ad un problema parabolico con frontiera libera. 


Le opzioni di tipo asiatico si differenziano da quelle di tipo europeo per il 
fatto che il loro valore al tempo t = T dipende non solamente da S (T), ma 
anche dalla sua media calcolata sull’intervallo [0,7], come risulta dai seguenti 
esempi: 


V(T,S(T)) = max SC7) - z d SIAdr0} ; (2.1) 
V(T, SO) = mex {i | S(r)dr — E, o} ; (2.2) 


V(T, S(T)) = max {st) isp È i log(S(r))dr) 0} . (2.3) 


V(T,S(T)) = max {esp È J i log(S(7))dr) Li 0}. (2.4) 
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Osserviamo che i processi stocastici che descrivono il valore di queste opzioni 
presentano caratteristiche non markoviane, è pertanto opportuno introdurre una 
nuova variabile A che tenga conto della storia del sitema e considerare V come 
funzione delle tre variabili indipendenti (t,5, A). Negli esempi (2.1) e (2.2) 
appena visti è naturale scegliere 


t 
A(t) = f S(7)dr, 
0 
negli esempi (2.3) e (2.4) invece 


A(1) = J log(S(7))dr. 


Procedendo sostanzialmente come nel caso delle opzioni di tipo europeo si ot- 
tiene che V deve verificare la seguente equazione differenziale 


1 30202V OVLOV gIV 


nel caso della media aritmetica, 


1,20 2V | g0V 9V IV 
37°S D52 +57 + DI +rSag TV = 0. (2.6) 


nel caso della media geometrica. 

Il prezzo dell’opzione risulta pertanto determinato dall’equazione differen- 
ziale (2.5) (oppure (2.6)) e da una delle condizioni finali precedentemente elen- 
cate. 


Solamente in alcuni casi particolari (più precisamente quando si considera il 
dato finale (2.3)) è nota un'espressione esplicita per la soluzione del problema 
sopra descritto; nel caso della media aritmetica, che corrisponde alle opzioni 
effettivamente trattate nei mercati, le cose diventano considerevolmente più dif- 
ficili. Ricordo che, sotto opportune condizioni, si è cercato di approssimare la 
densità di probabilità del processo stocastico corrispondente all’equazione (2.5) 
con la densità corrispondente alla (2.6), controllando l’errore per mezzo della 
conoscenza dei suoi momenti. 

Quando si considera direttamente l'equazione (2.5) non si conosce un’espres- 
sione esplicita per la soluzione del corrispondente problema di Cauchy e si è co- 
stretti a ricorrere a metodi numerici, che possono mirare a simulare il processo 
stocastico (come ad esempio nei lavori di Kemna e Vorst [8], di Ruttiens [15] e 
di Caverhill e Clewlow [5]), oppure possono essere basati sull’approssimazione 
numerica della soluzione del problema di Cauchy (si vedano i lavori di Barra- 
quand e Pudet [3], Barles, Daher e Romano [2], Wilmott, Dewynne ed Howison 
[17], Dewynne e Wilmott [6], Rogers e Shi [14]). Nel nostro lavoro abbiamo 
scelto quest’ultima strada, cercando dapprima di ottenere risultati di esistenza, 
unicità e regolarità per le soluzioni del problema di Cauchy, per poi costruire 
metodi di approssimazione numerica alle differenze finite. 
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3 Media geometrica 


Consideriamo dapprima il caso della media geometrica. 
Se poniamo x = Y2 log(S), y= YZ A e 
2r-0? 2r a? 2 sa (dì 
u(2, y;t) = e E787+ (7) j V (7 Lat t,ev?, cu) ) 
vV2 

allora l’equazione (2.6) risulta equivalente all’equazione considerata da Kolmo- 
gorov in [16, pag.167] 

du du _ du 

dx? Oy dt : Vili 
che costituisce il prototipo della classe di operatori studiata in [9] dall’autore 
e da Lanconelli, mentre le condizioni finali (2.3) e (2.4) sì trasformano nelle 
condizioni iniziali 


u(7,y,0) = e Ere” max {eF - eFF,0} (3.2) 
e 
ar-o? CA 
u(2,y,0)=e avo * max fer _ E,0}. (3.3) 


Ricordiamo alcuni risultati noti per il problema di Cauchy 


(3.4) 


Sur =%, per (2,y,t) ERx Rx]O, TT, 
u(2,y,0) = (x,y) per (1, y) E R?. 


Innanzitutto la funzione 


v3 - è? 3 t- ? 
rt = pie E (240) 
(3.5) 


per t > 7, I(2,y,t;€, r)=0 pert<r,è una soluzione fondamentale di [8.3], 
quindi, se 9 è una funzione continua su R? e tale che 


J e7(°+8°)p(2, 4) dr dy< 0%, 
R2 


per una costante arbitraria c, allora la funzione definita da 
u(2,y,t) = JFEnten06l6n dé dn (3.6) 


è soluzione del problema di Cauchy (3.4). Vale inoltre il seguente risultato di 
unicità: se u e v sono due soluzioni di (3.4) € verificano la seguente condizione 


J e-0E"+9°) |u(x, 9,4) 2 v(x2,y,4)| dx dy di < 00, 


ST 
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per una costante arbitraria c, allora u = v (sì veda [12]). 


I risultati sopra enunciati si possono estendere al caso in cui i coefficienti 
o ed r sono funzioni limitate e localmente hélderiane della variabile t, con o 
strettamente positiva. Se infatti R indica una primitiva di r e se poniamo 


u(z,y,t)= eRO vr te”,U): 


allora l’equazione differenziale (2.6) risulta equivalente alla 


o(t)? du o(t)?\ du du _ du - 
2 dr? + (r0 son Do) —- (3.7) 


de dI di 
che rientra tra quelle studiate in [11]. Naturalmente nella formula (3.6) si deve 
sostituire la funzione F con la soluzione fondamentale dell’equazione (3.7). 


Per quanto riguarda l’approssimazione numerica della soluzione di (3.4) pos- 
siamo adattare alcuni risultati di [13], dove è stato studiato un problema di valori 
al contorno per un aperto limitato di R3. Definiamo l’operatore “discreto” 


Lou(z,g,t)= (ut -uley+ér,t-8), 


u(e — h,y+òx,t-6d) — 2u(2,y+dr,t-J) +ulr+h,y+d6x,t-6ò) 
TTT fr e I ’ 
che approssima l’operatore L nel senso che soddisfa alla seguente condizione 


LGu= Lu+(h°+6)O(h,6) 


(qui O(A, 6) indica una funzione limitata). L’operatore LG risulta ben definito 
sulla griglia 


G={(jAx,kKAy,nA,}) €R?:j,k,n e z}, 
con Az=h,A.=lée Ay= hd e vale il seguente Teorema di approssimazione: 


TEOREMA 3.1 Sia u soluzione di (3.4), ug soluzione di Lgue = 0 inGn 
ST,UG=  inGN {t = 0}. AWora, per ogni €1,€2 > 0 e per ogni compatto H 
di Sr esiste una griglia G, verificante la seguente condizione di stabilità 


Ai 1 
gii 
(Ax)? 72° 
tale che 
” n i) — t){< A 
(retina (3:91) = vole, yi) <a 
u uG(t+Ax,y,t)- ug(2, Yy,t) 
Cda yi) RETI e). 
(ev) EnH dx (2,9,8) Az nia 


Risulta inoltre e; = O((Ax)?) ed e2= O(Ax). 
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Osserviamo che il metodo fornisce anche un’approssimazione della quantità 
A definita in (1.5) che nel caso delle opzioni di tipo europeo viene effettivamente 
utilizzata dagli operatori nella gestione del portafoglio che era stato indicato con 
II. 

Ricordiamo infine che, per utilizzare gli stessi metodi di approssimazione 
numerica nel caso dei coefficienti o ed r variabili, è opportuno richiedere che 
questi ultimi siano derivabili ed abbiano derivata localmente hòlderiana. 


4. Media aritmetica 


Prima di studiare il problema di Cauchy corrispondente al caso della media 
aritmetica, sserviamo che il cambiamento di variabile r = log(.5) che abbiamo 
effettuato nel paragrafo precedente trasforma l’equazione differenziale di (2.5) 
in un’equazione del tipo seguente 


du du _ du 


ate du di 


"i (4.1) 


che presenta delle difficoltà decisamente diverse rispetto alla (3.1). Sappiamo che 
l’equazione (4.1) conserva molte proprietà di tipo “locale” dell’equazione di Kol- 
mogorov (3.1): ad esempio rientra anch’essa nella classe studiata da Hormander 
in [7], e quindi ogni sua soluzione è necessariamente una funzione di classe C'°°. 
Proprietà di tipo “globale” relative alle soluzioni di (4.1) sono molto più diffi- 
cili da ottenere, principalmente per il fatto che il coefficiente della derivata qu 
cresce esponenzialmente. Per fare riferimento in maniera più diretta a risultati 
noti, ci è sembrato preferibile lavorare sull’equazione nella forma.con cui appare 
nella (2.5). 


Se poniamo 


allora l’equazione differenziale (2.5) risulta equivalente alla 


20° du _ du, 


° 92? * °9y De: 
di conseguenza il problema di Cauchy (2.5) con dato (2.1) (o (2.2)) diventa 
2°94 + 90 = Qu per (r,y,t) € Qr, (42) 
u(2,v,0) = (x,y) per (2,9) ER+TxR, 


dove abbiamo posto Q7 = RT x Rx]0, 2T[ e 


g(xc,y) = e" max (è - Aol 0) (4.3) 
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nel caso di (2.1) (siamo interessati a valori positivi di y), 


do E,0) (4.4) 


., —_ mn Xx 
P(2,y) = e” ma ($ 


nel caso di (2.2). 
Enunciamo ora un risultato di esistenza di soluzioni per il problema (4.2). 
A tal fine faremo uso della nozione di sopra e sottosoluzione. Diremo che u è 
soprasoluzione del problema (4.2) se 
2204 4 p@ù < du per (x,y,t) € Qr, (4.5) 
u(x,Y,0) > p(2,y) per (2,9) € R+xR 


e che u è sottosoluzione di (4.2) se vale la condizione (4.5) con le due disugua- 
glianze invertite. 


TEOREMA 4.1 Siano g € C(RtxR),u u rispettivamente soprasoluzione e sot- 
tosoluzione del problema (4.2) e supponiamo che u < © in Qr. Allora esiste 
una soluzione u di (4.2) tale che 


u<u<U inQr. (4.6) 


La prova del Teorema sopra enunciato si basa sui risultati di esistenza provati 
da Manfredini in [10] per il problema di Dirichlet su un aperto limitato e su 
un procedimento diagonale. Osserviamo che per ottenere una soluzione del 
problema relativo al dato iniziale (4.3), è sufficiente scegliere u = 0 ed u(2,y,t) = 
x"+1e, mentre per il problema relativo al dato (4.4), u= 0 ed T(2,y,t) = 
debtam [72 +y2,conk=m?+2m+2. 

In entrambi i casi risulta u(0,y,t) = (0, y,t) = 0, quindi abbiamo implici- 
tamente richiesto che u(0, y,t) = 0 per ogni y, t; di conseguenza, per il Teorema 
B di Aronson e Besala [1], u è l’unica soluzione che verifica (4.6). 

Nel caso particolare del problema (4.2), l’ipotesi u(0,y,t) = 0 non è effetti- 
vamente necessaria, come mostrato nel seguente Teorema. Definiamo prelimi- 
narmente su R+ x RR la seguente funzione 


(2,9) = log (2° +y° +1) — log(x) (4.7) 

ed osserviamo che risulta 
V(a,y) co, (x,y) + 0, 
Il(2,y)ll +00 (2,4) +(0,Y0) 
per ogni yo E R. 
TEOREMA 4.2 Sia u € C(Qr), tale che Lu < 0 in Qr, che u(x,y,0) > 0 per 
ogni (x,y) E Rtx R e supponiamo che esistano due costanti positive M e k tali 
che risulti , 
u(2,y,t) > -Mebt (9), 

per ogni (x,y,t) € Q7 (4 indica la funzione definita in (4.7)). Allora u> 0 in 
Or. 
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Per quanto riguarda i metodi di approssimazione numerica i nostri risul- 
tati non sono ancora completi. Osserviamo che la soluzione del problema (4.2) 
viene ottenuta come limite di una successione di soluzioni di problemi relativi 
ad aperti limitati e quindi è sufficiente approssimare una di queste soluzioni. Di 
conseguenza è possibile trascurare il comportamento della soluzione “all’infini- 
to” e si possono adattare senza difficoltà i metodi numerici proposti nella nota 
[13]. Certamente questi risultati non sono ancora molto soddisfacenti, in quanto 
non forniscono una stima esplicita dell’errore, ma almeno indicano che queste 
tecniche possono effettivamente fornire una buona approssimazione numerica. 


Per concludere, osserviamo che la seguente equazione differenziale stocastica 
dS = oS°dW + uSdt, 


con 0 < 0 < 1, è stata proposta come modello più realistico per la descrizione 
dei mercati. In questo contesto l’equazione differenziale che determina il prezzo 
dell’opzione diventa 


dV + giV + 4 ba 
os? dA dt Chi 

le tecniche indicate in questo paragrafo si adattano senza alcuna difficoltà a 
questo tipo di problemi, anche nel caso di coefficienti 0, p ed r dipendenti da 1. 


50?s% rV=0; 
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